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CHAPTER  I 
INTRODUCTION 

In  recent  years,  the  transmission  voltage  used  in 
electrical  power  systems  has  climbed  to  700  KV  and  some  con- 
siderations have  been  given  to  using  1,000  KV.   In  such  EHV 
(extra  high  voltage)  systems,  high  performance  air  blast 
circuit  breakers  are  being  used  as  a  vital  element  for  pro- 
tection of  the  power  'transmission  system. 

This  type  of  circuit  breaker  has  economic  and  other 
advantages  over  the  bulk  oil  type  circuit  breakers  which  have 
been  utilized  for  many  years  world  wide. 

In  an  air  blast  circuit  breaker,  compressed  air  stored 
in  a  fixed  volume  supply  tank  is  the  main  operating  vital. 
The  compressed  air  serves  three  purposes: 

1.  It  provides  the  energy  for  opening  or  closing  the 
circuit  breaker  contacts  on  command  from  a  control 
mechanism. 

2.  It  blasts  away  the  products  of  the  arc  and  cools  the 
arc  column  during  breaker  opening  (interruption). 

3.  It  provides  a  dielectric  medium  across  the  circuit 
breaker  contacts,  when  they  are  open,  and  also 
between  live  parts  and  ground. 

For  a  circuit  breaker,  fast  and  positive  interrupting 
operation  is  considered  to  be  of  utmost  importance.   Usually, 
it  demands  that  interruption  must  occur  within  several 


hundredths  of  a  second  in  case  of  a  short  circuit  in  the 
transmission  system  due  to  faults  such  as  falling  trees, 

lightning,  etc.   Comprehensive  discussion  and  study  in  recent 

p 

years  has  shown  a  desire  for  breakers  with  2  cycles  (e.g.  ^q 

seconds)  interrupting  time  rather  than  the  present  American 
standard  of  3  cycles  (1)*.   Benefits  of  the  shorter  inter- 
rupting time  includes  reduced  arcing  damage  to  lines  and 
cables  from  faults  and  higher  line  loading  capability  through 
improved  stability  because  of  shorter  fault  duration  time. 

In  designing  a  pneumatic  control  system  for  an  air  blast 
circuit  breaker,  a  primary  consideration  for  its  mechanical 
operation  and  interrupting  time  is  the  operating  pressure 
to  be  used.   Contact  velocities  and  valve  operation  are  mainly 
a  function  of  the  pressure  of  the  air.   Beyond  the  consider- 
ation of  pressure,  the  only  requirement  is  that  the  air  be 
relatively  clean  and  free  from  liquid  water  and  solid  foreign 
objects.   This  requirement  is  imperative  for  practical, 
efficient  and  dependable  operation  of  an  air  blast  circuit 
breaker. 


♦Numbers  in  the  parenthese  refer  to  items  of  references. 


CHAPTER  II 

DESCRIPTION  AND  ANALYSIS  OF  A  PNEUMATIC  CONTROL 
SYSTEM  FOR  AN  AIR  BUST  CIRCUIT  BREAKER 


In  this  report,  interests  are  focused  upon  the  response 

of  a  piston  type  actuator  driving  a  spring  loaded  set  of 

i 
circuit  breaker  contacts. 

The  complete  system  is  shown  schematically  in  Fig.  1 
and  consists  of  a  fixed  volume  supply  tank,  a  piston  type 
actuator,  and  an  on-off  solenoid  type  control  valve. 

It  is  assumed  that  the  actuating  system  controls  the 
opening  and  closing  of  a  set  of  circuit  breaker  contacts 
directly,  so  that  considerations  of  linkage  delays  and 
backlash  can  be  eliminated.   In  an  actual  circuit  breaker, 
the  number  of  sets  of  contact  is  determined  by  the  current 
capacity  requirement  of  the  power  transmission  system.   For 
instance,  in  an  Allis-Chalmers  Type  ABM  air  blast  circuit 
breaker  (2),  there  may  be  as  many  as  six  sets  of  similar 
contacts  in  each  phase  of  the  breaker. 

As  described  in  Chapter  I  of  this  report,  in  an  actual 
circuit  breaker,  the  compressed  air  stored  in  a  supply  tank 
serves  three  functions  during  breaker  operation.   Therefore, 
the  air  pressure  in  the  supply  tank  undergoes  a  considerable 
drop  during  an  actuating  cycle.   For  example,  when  a  type 
ABM-500-38,000  circuit  breaker  executes  an  opening  (inter- 
ruption) cycle,  a  supply  tank  pressure  drop  of  approximately 
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Fig.  1.   Schematic  Representation  of  an  Actuating 
System  for  an  Air  Blast  Circuit  Breaker. 


36  psi  occurs.   In  this  report,  the  effect  of  changing  supply 
tank  pressure  on  actuator  response  time  will  be  investigated. 

In  order  to  introduce  different  drop  in  supply  tank 
pressure  during  opening  of'  the  breaker  contacts,  different 
tank  volumes  are  assumed. 

The  following  parameter  values  for  system  components 
were  selected  for  study  in  this  report.   These  values  were 
selected  as  being  typical.  However,  the  analysis  techniques 
developed  in  this  report  can  be  applied  in  general  for  any 
desired  combination  of  values  of  system  parameters. 

The  parameter  values  selected  were: 

1.  Vt  -  1,000  in3,  500  in3,  200  in3  and  70  in3. 

2.  A0  =  0.194  in2,  0.097  in2,  0.04^5  in2,  0.0277  in2 

and  0.0194  in2. 

3.  Kd  -  0.95. 

4.  A  =  3-14  in2. 

5.  Lc  -  8  in. 

6.  Cn  =  0.5  in. 

o  • 

7.  b  =  0.75  in. 
3.  IVL  »  1.126  lbm. 
9.   ML  -  10  lbm. 

10.  B =  1.5  lbf-sec/in. 

11.  Ks  =  291  lbf/in. 

12.  Cr  =  3  in.,  3.5  in.,  4  in. 
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Following  assumptions  were  used  in  the  analysis  of 

the  system: 

1.  The  working  fluid,  air,  is  a  perfect  gas  (4). 

K  =  1.4 

R  =  53  .,34  ft-lbf/lbm-°R 

2.  The  working  fluid  has  constant  specific  heats. 

3.  The  kinetic  energy  of  the  working  fluid  is 
negligible. 

4.  Flow  velocity  is  high  enough  to  consider  flow 
processes  as  adiabatic. 

5.  Coulomb  type  friction  is  negligible. 

6.  The  supply  tank  actuating  cylinder,  control  valve 
and  interconnecting  lines  are  rigid. 

7.  The  discharge  coefficient  of  the  control  valve  is 
constant . 

8.  Pressure  wave  effects  in  the  working  fluid  are 
negligible. 

9.  At  time  t  =  0,  the  solenoid  valve  is  closed;  at  time 
t  =  0+,  the  solenoid  valve  is  wide  open. 

Initial  conditions  assumed  were: 

1.  t  =  0+ 

2.  Pt(0)  -  *370  psig.  =  384.7  psia. 

3.  Pc(0)  =  Patm  =14.7  psia. 


*This  is  the  supply  tank  pressure  used  in  Allis-Chalmers 
type  ABM  air  blast  circuit  breakers  (2). 


4.  et(o)  -  530  °R. 

5.  6C(0)  -  530  °R. 

6.  ft.(0)  =  0.001136  lbm/in3. 

7.  p  (0)  -  0.0000434  lbm/in3. 

8.  Pt(0)  -  0  lbf/in2-sec. 

9.  Pc(0)  -  0  lbf/in2-sec. 

10.  C(0)  -  0  in. 

11.  C(0)  -  0  in/sec. 

12.  *C(0)  =  in/sec2. 

A.  Flow  Through  the  Solenoid  Control  Valve 

Consider  opening  of  the  breaker  contacts.  At  time 
t  -  0+,  the  solenoid  valve  is  wide  open,  and  the  pressure 
ratio  is: 

^__  .  j-W  =  0.0382 
Pt(0)   384.7 

This  ratio,  0.0382,  is  less  than  0.528,  the  criterion  for 
distinguishing  between  critical  and  subcritical  flow  (3)  (4) 
(5).  Therefore,  the  flow  through  the  control  valve  at  the 
beginning  of  an  opening  cycle  is  critical. 

The  equation  for  critical  flow  through  the  solenoid 
valve  is: 


2  ,  K+l 


VT(t)  -  KdA   /  Kljfj)  ^  p;<t)  Pt(t)  (2.1) 
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In  equation  (2.1),  W'(t)  has  the  dimensions  lbf sec/in. 
P£(t)  has  the  dimensions  lbf  'sec2/iiA.   Converting  the 
dimension  of  W(t)  to  lbm/sec,  the  dimensions  of  Pt(t)  to 
lbm/in^,  and  substituting  K  =  1.4,  equation  (2.1)  becomes: 


W(t)  =  KdAQ   /180.  5   Pt(t)  Pt(t) 


(2.2) 


Equation  (2.2)  is  applicable  until  the  following  condition 
prevails: 

Pc(t) 

p~(7)  >  0.528 


Then  the  flow  through  the  solenoid  control  valve  becomes 
subcritical.   The  subcritical  flow  equation  is: 


Pjt) 


^t)  =  KdAo    W;(t)  V^to 


K-Hr 
K 


Pc(t) 


-2/K 


(2.3) 
Converting  the  dimensions  of  W'(t)  and  ^(t)  into 
lbm/sec  and  lbm/in^  respectively,  and  substituting  K  =  1.4, 
equation  (2.3)  becomes: 


W(t)  =  51.9KdA0     Pt(t)Pt(t) 


r  pp,(t) 


1 


1 0.286' 


Pt.(t) 


•Pc(t) 


PTTtT 


1.428 


(2.4) 


B.   Rate  of  Mass  Flow  from  the  Supply  Tank 

If  the  supply  tank  is  considered  as  a  control  volume, 
the  equation  which  defines  the  rate  of  mass  flow  from  the 
tank  is: 

«»-M»i;ifl(,) 

Substituting  and  rearranging  yields: 
-1.4Pt(t)V/(t) 

?^]-     vtet(t)       ■  (2-6) 

C.   Rate  of  Mass  Flow  into  the  Actuating  Cylinder 

Considering  the  actuating  cylinder  as  a  control 
volume,  the  rate  of  mass  flow  into  the  cylinder  is: 

W(t)  =  rAC0+AC(t)l  ?cU)  e°it]   +  Pc(t)AC(t)       (2.7) 
L  °      J  KPc(t) 

Substituting  K  =  1.4  and  rearranging  gives: 

* 

.       1.4[w(t)-  pc(t)AC(t)J  Pc(t) 

P  (t)  =  != LZ. -J 2 (2.8) 


A  Pc(t)  C0+C(t)J 


D.  Forces  Acting  on  the  Actuator  Piston 

Summation  of  the  forces  acting  on  the  piston,  yields 
the  force  balance  equation: 

Pc(t)A-PatmA-(ML+Ml)C(t)-BpC(t)-KsC(t)  =  0        (2.9) 
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The  dimension  of  M£  and  Ml  in  equation  (2.9)  is  lbf. 
sec  /in.   Converting  the  dimensions  of  Mc  and  Mp  to  lbm, 
then  rearranging,  equation  (2.9)  becomes: 

386{[pc(t)-Patm]A-BpC(t)-KsC(t)) 
C(t)  -  M^  u.iu; 

It  is  noticed  that  the  above  derived  equations  (2.1) 
through  (2.9)  are  all  nonlinear  except  equation  (2.10).   Ob- 
viously the  solution  of  these  equations  needs  special  treatment 

E.  Approximate  Solution  to  Obtain  the  Response 
of  the  Nonlinear  System  (6)  (7) 

If  the  values  of  P(t),  Pc(t),  C(t),*C(t)  are  considered 
to  be  constant  for  a  small  time  interval  At  which  is  much 
smaller  than  the  dominant  response  time  of  the  system,  the 
following  equations  apply  for  calculating  the  values  of 
Pt,  Pc,  t   and  C  at  t  =  0+At : 

Pt(0+At)  =  Pt(0)  +  Pt(0)At  .    (2.11) 

Pc(0+At)  =  Pc(0)  +  Pc(0)At  (2.12) 

C(0+At)  =  C(0)  +*C*(0)At  (2.13) 

C(0+At)  -  C(0)  +  C(0)At  (2.14) 

When  a  perfect  gas  undergoes  an  adiabatic  expansion,  such 
as  occurs  in  the  solenoid  valve,  the  following  equations  apply: 


. 
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•Pt(0+At)  =  Pt(0) 


pc(o+At)  =  ec(o) 


pt(o+^t) 

Lpt(o)  : 


0.714 


Pr(0+At) 


-10.714 


.P) 


(2.15) 


(2.16) 


Values  calculated  using  equations  (2.11)  thru  (2.16)  can  be 
substituted  into  equations  (2.2),  (2.6),  (2. 8)  and  (2.10)  to 
compute  W(0+At),  Pt(0+At),  Pc(0+At)  and*C*(0+At)  respectively. 

Starting  at  t  =  0+2At  and  continuing  in  increments  of 
At,  the  Runge-Kutta  approximation  can  be  utilized  in  order 
to  obtain  better  accuracy  in  this  type  of  approximate  solu- 
tion. 

Using  the  Runge-Kutta  method  all  equations  can  be 
generalized  for  t  =  0+nAt : 

Pt(0+n*t)  =  pl  [n+(n-n  W.]  vfet  (0+"*t  )+Pt  [o-Kn-D^At      {2>1?) 

p  (*»*)  -  p.  ro+(^iH+fc((Hnot)+*°  ^'"-i'A4At  (2.18) 

c  c  L  J  2     ■ 

.  r                   -,  {c(0+nAt)+cT6+(n-l)Atl"Ut 
C(0+nAt)   =  C  |0+{n-l)AtJ+i -*= ^ (2.19) 

P     ,        %     i   {c(0+nAt)+C  fb+(n-l)AtlUt 
)   -  C  [0+(n-l)AtJ+i — 1= ii (2.20) 


C(0+nAt 


Pt(0+nAt)   =  p     0+(n-l)At 


Pt(0+nAt) 
Pt  [0+(n-l)At]_ 


0.714 


P  (0+nAt)   =  Cc  0+(n-l)At  < 


Pc(0+nAt) 


Pc  jO+(n-l)At]J 


0.714 


(2.21) 


(2.22) 
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W( 0+nAt)  =  KdA0  /  180.5  Pt  (0+n^t )?t (0+nAt ) 
(If  critical  flow) 


(2.23) 


W(0+nAt)=51.9KdAo  /  Pt  (0+nAt)Pt  (0+nAt ) <h-  p^(Q+nAt) 


Pc (0+nAt ) 


0.236 


Pr  (0+nAt) 


Pt  (0+nAt) 


1.428  ' 


(2.24) 


(If  subcritical  flow) 


a    /„      ,    i       -1.4Pt (0+nAt )W  0+nAt 

Pt   0+nAt     -  V  . -J - 

t  VtPt  0+nAt 


(2.25) 


•    ,_          ,      1.4  (w ( 0+nAt  )-Pr  (0+nAt  )AC  (0+nAt  )~|pp(0+nAt) 
Pc(°+ntt)-  — i Aec(0^t)[Co(0+nAt)j  '2'26> 


C(  0+nAt  )  = 


386-![pc  (0+nAt  )-Patm  ]  A-BpC  (0+nAt  )-KsC  (0+nAt  )~j> 


VL  +   MT 


% 


(2.27) 


Based  upon  this  set  of  equations,  a  flow  diagram  for 
solving  the  equations  iteratively  using  an  IBM  360  Computer 
was  determined  (see  appendix  A).  Then  a  FORTRAN  program  was 
written  (see  appendix  B)  for  the  computation  of  Pt,.Pc,  C 
and  C  with  respect  to  time  t  as  a  variable. 
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CHAPTER  III 
DISCUSSION  OF  RESULTS 
A.  Actuator  Response  Characteristics 

Initially,  data  were  computed  for  Vt  =  1000  irK, 
500  in^  and  200  in3  respectively.   (Fig.  2,  3  and  4).   Cal- 
culated data  included  piston  displacement,  piston  velocity, 
supply  tank  pressure  and  cylinder  pressure  as  functions  of 
time  following  a  step  change  in  solenoid  control  valve 
area  from  0  to  AQ  at  time  t  =  0+.  _ 

Fig.  2  shows  that,  for  a  supply  tank  volume  of  1000  in->, 
a  solenoid  valve  metering  area  of  0.0194  in^  and  an  initial 
supply  tank  pressure  of  384*7  psia,  the  supply  tank  pressure 
dropped  to  3^1.6  psia,  after  an  elapsed  time  of  0.072  sec. 
The  total  supply  tank  pressure  drop  for  the  complete 
actuating  cycle  is  only  3.1  psi. 

Mote  the  cylinder  pressure  curve,  Pc (t ) ,  and  the  piston 
velocity  curve,  C(t).  At  the  beginning,  both  Pc(t)  and  • 
C(t)  increase  very  rapidly.   This  tendency  persists  until 
t  =  0.0055  sec.  At  this  instant,  the  volume  in  the  cylinder 
is  increasing  very  rapidly,  approximately  25.13  in^/sec. 
Apparently  because  of  the  rather  small  value  of  A0,  the  air 
flow  rate  entering  the  cylinder  is  not  sufficient  to  keep 
cylinder  pressure  increasing  and  it  starts  to  drop.   The 
decrease  in  Pc(t)  causes  C(t)  to  slow  down.   C(t)  has  its 
peak  value  at  t  =  0.0135  sec.  After  this  C(t)  keeps 
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decreasing  and  the  interaction  causes  Pc(t)  to  increase 
until  it  equals  the  tank  pressure  Pt(t)  at  the  time  that  the 
piston  comes  to  rest. 

Consider  the  piston  displacement  curve,  C(t),  it  in- 
creases very  rapidly  during  the  early  stages.   It  increases 
approximately  linearly  between  t  =  0.01  sec  and  t  =  0.02 
sec.  At  t  =  0.072  sec  the  system  has  very  nearly  reached 
its  steady  state  condition  and  C(0.072)  =  4.05  in. 

The  accuracy  of  the  calculated  results  can  be  checked 
as  follows: 

At  t  =  0.072  sec  C^'O  and  c'^0,  therefore,,  the  force 
balance  equation  for  the  piston  is: 

[Pc(0.072)-Patm]  A-KSC(0.072)  =  0  (3-D 

or 

P  (0.072)-PQa.„  ,  •' 

C(0.072)  -  -c atm  (3.2) 

Ks 

Substituting  ?Q (0.072)  -  3^1.6  psia,  Patm  =  14-7  psia, 

A  =  3.14  in2  and  Ks  =  291  lbf/in. ,  into  equation  (3.2)  gives: 

C(0.072)  -  3.9S  in. 

Compared  to  the  approximate  result  of  the  numerical  solution 
the  relative  error  is: 

This  error  is  tolerable  in  this  application. 
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For  a  tank  volume  of  500  in^ ,  a  valve  metering  area  of 
0.0194  in2  and  an  initial  supply  tank  pressure  of  3&4.7  psia, 
Figure  3  shows  that  the  supply  tank  pressure  dropped  384.7- 
37S.7  -  6  psi. 

For  Vt  =  200  in3,  A0  =  0.0194  in2  and  Pt(0)  =  3^4-7 
psia,  Figure  4  shows  that  the  supply  tank  pressure  drop  for 
a  complete  actuating  cycle  is  3&4. 7-370.6  =  14. 1  psi. 

It  can  be  observed  that  the  Pc(t),  C(t),  and  C(t) 
curves  for  both  Fig.  3  and  Fig.  4  are  very  similar  to  the 
same  respective  curves  in  Fig.  2. 


B.  Actuator  Displacement  Characteristics 


As  discussed  in  Chapter  I  of  this  report,  for  an  air 
blast  circuit  breaker,  the  interrupting  time'  (the  time  re- 
quired for  opening  of  the  contacts)  is  desired  to  be  less 
than  0.033  sec.  (2  cycles).  In  other  words,  the  time  re- 
quired for  the  piston  (directly  connected  to  the  breaker 
contacts)  to  move  a  certain  distance  is  of  prime  interest 
in  this  report. 

To  emphasize  this  interest,  a  distance  quantity  Cr  is 
hereby  defined.   Cr  is  defined  as  the  necessary  piston  dis- 
placement so  as  to  acquire  complete  opening  of  the  contacts. 
For  the  purposes  of  this  report,  Cr  will  be  assumed  to  be 
3,  3-5,  or  5  inches.   Next,  a  time  quantity  tr  is  defined  as 
the  required  response  time  for  the  piston  to  move  Cr  =  3, 
3.5  or  4  inches.   These  2  quantities  will  be  used  as  a  basis 
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for  investigating  and  comparing  the  approximate  actuator  dis- 
placement vs  time  curves  calculated  when  certain  parameters 
of  the  actuating  system  are  varied. 

In  order  to  obtain  a  response  characteristic  for  a 
simulated  tank  pressure  drop  of  36  psi  (which  corresponds 
to  the  supply  tank  pressure  drop  experienced  in  a  complete 
opening  of  an  Allis-Chalmers  type  ABM  air  blast  circuit 
breaker),  it  was  necessary  to  select  a  supply  tank  volume 
of  70  in^.  Actuator  displacement  vs  time  curves  for  this  case 
are  shown  in  Fig.  5  for  A0  =  0.194  in  ,  and  in  Fig.  6  for 
A0  =  0.019A  in2.   For  convenience  of  comparison,  the  actuator 
displacement  curve  for  a  tank  volume  of  1000  in^  is  also 
plotted  on  each  figure. 

Since  piston  displacement  up  to  4  inches  is  of  most 
interest,  the  displacement  curves  in  Fig.  5  (as  well  as  in 
Fig.  6,  7,  8,  etc.)  are  not  completely  presented  for  the 
purpose  of  simplicity. 

In  Fig.  5,  for  a  control  valve  area  of  AQ  =  0.0194  in2 
and  a  supply  tank  volume  of  V*t  =  1000  in^,  the  pressure  drop 
during  the  complete  actuating  cycle  is  only  3.1  psi.,  as  was 
previously  noted. 

In  the  same  figure,  for  the  same  control  valve  area  and 
a  supply  tank  volume  of  70  in^ .   The  pressure  drop  during  the 
complete  actuating  cycle  is  now  35.8  psi.   This  is  approximately 
the  same  amount  of  pressure  drop  which  occurs  in  the  Allis- 
Chalmer  type  ABM  circuit  breaker  during  an  operating  cycle. 
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Fig.  5.   Actuator  Displacement  Characteristics 
A0  =  0.0194  in2. 
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Fig.  6.  Actuator  Displacement  Characteristics 
A0  -  0.194  in2. 
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Comparing  the  response  time  tr  for  a  displacement  of 
C  -  3  in.  the  time  difference  between  these  two  pressure 
drops  is  only  0.0005  sec. 

For  Cr  =  3.5  in.,  the  difference  between  the  two  response 

times  is  0.0025  sec. 

It  can  be  observed  that  as  time  increases  the  difference 
in  time  required  to  reach  a  given  displacement  becomes  more 
significant. 

In  Fig.  6  for  a  solenoid  control  valve  area  AQ  -   0.194 
in2,  the  respective  pressure  drop 'of  the  1000  in3  tank  and 
70  in^  tank  are  the  same  as  that  of  Fig.  5.   The  response  time 
difference  for  Cr  =  3  in.  is  also  0.0005  sec. 

For  Cr  »=  3.5  in.,  the  tr  difference  is  0.0006  sec  which 
is  smaller  than  that  of  Fig.  5-  For  Cr  -  4  in.,  the  tr  dif- 
ference becomes  appreciable,  i.e.,  tr  ■  0.001  sec. 

In  summary,  it  can  be  said  that  for  supply  tank  pressure 
drops  as  large  as  that  experienced  in  the  Allis-Chalmers 
type  ABM  circuit  breaker,  during  its  opening  operation 
cycle,  there  is  no  significant  difference  in  response  time 
for  actuator  displacements  up  to  3.0  inches  (for  the  assumed 
system) . 

C.   Effect  of  Changing  Control  Valve  Metering 
Area  on  Response  Characteristics 

In  Fig.  7  five  piston  displacement  curves  are  shown  for 
a  system  with  a  supply  tank  volume  of  1000  in^  and  solenoid 
control  valve  areas  of  0.194  (A0  =  1),  0.097  (A"0  =  0.5), 
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Fig.  7.     Actuator  Displacement  Characteristics 


Vt.  =  1000  in3. 
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0.04^5  (A0  =  0.25),  0.0277  (A0  =  0.143)  and  0.0194  (AQ  =  0.1) 

2 
in  respectively.    i 

Observe  the  curves  plotted  in  Fig.  7  for  A0  -  1  the 
response  time  for  C  ■  3  in.  is  0.0155  sec.   If  the  solenoid 
valve  area  is  reduced  by  one  half,  AQ  ■  0.5,  tr  is  increased 
slightly  to  0.016  sec. 

If  the  valve  area  is  reduced  to  A0  =  0.143,  tr  is  in- 
creased to  0.0195  sec.   Finally,  if  the  valve  area  is  reduced 
to  A0  -  0.1,  tr  increases  to  0.026  sec. 

It  is  apparent  that  the  area  A0  of  the  on-off  control 
valve  is  a  significant  parameter  as  far  as  response  time  is 
concerned.  The  larger  the  solenoid  valve  area,  the  faster 
the  piston  moves,  and  thereby  the  shorter  the  response  time 

V 

Therefore,  the  response  time  t  of  the  actuator  can  be 
increased  or  decreased  depending  on  the  area  selected  for 
the  on-off  control  valve.   But  one  thing  must  be  pointed  out 
here,  that  a  limit  for  increasing  A0  exists.  As  the  area 
A0  (Fig.  7)  is  increased  from  0.25  to  1.0.   The  response 
time  tr  is  only  reduced  from  0.0192  sec.  to  0.0155  sec, 
(0.0037  sec),  while  increasing  A0  from  0.1  to  0.25,  reduced 
tr  from  0.026  sec  to  0.0193  sec  (0.0093  sec). 

Similar  relationship  between  the  piston  displacement 
and  control  valve  area  can  be  observed  in  Fig.  &,   for  a 
supply  tank  volume  of  500  in^,  and  in  Fig.  9  for  a  supply 
tank  volume  of  200  in^ . 
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Fig.  8.  Actuator  Displacement  Characteristics 


Vt  =  500  in3. 
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Fig.  9..J  Actuator  Displacement  Characteristics 
Vt  =  200  in-3. 
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D.   Overshoot  Behavior  of  the  Piston  Displacement 

Fig.  10  shows  the  characteristic  response  curves  for  a 
supply  tank  volume  of  200  in^  and  a  solenoid  control  valve 
area  of  A0  =  0.5.   It  is  observed  that  both  the  displacement 
C(t)  curve  and  the  piston  velocity  C(t)  curve  overshoot 
considerably.   The  causes  .of  this  behavior  has  to  be  cleared. 

From  equation  (2.10)  in  Chapter  II.   The  following 
actuator  transfer  function  can  be  derived: 


C(s)        1 

F^T-  <M£+Mj)s2+BpS+Ks  (3'2) 

F(s)  is  the  Laplace  transform  of  the  pressure  force 
f(t)  acting  on  the  piston  of  the  actuator  and  is  defined: 

f(t)  «  [Pc(t)  -  PatJA  (3.3) 


The  damping  ratio  of  this  system  is 
Bd  1.5 


2  JKS(M.[+Mp)   2  /291(1Q+1-126) 


=  0.259  • 


In  Fig.  10  the  control  valve  metering  area  selected  is 
quite  large  and  it  can  be  observed  that  the  pressure  force 
Pc(t)  at  time  t  =  0+  is  very  nearly  constant  and  constant 
thereafter. 

It  is  approximately  a  step  input  for  an  underdamped 
second  order  system  having  a  damping  ratio  ^<  1.0.   The  dis- 
placement following  a  step  input  would  be  of  the  form 
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c(t)  = 


u) 


n 


u 


ni 


1-  co2 


■2 


r^uJnt     sinfcu^l-^    t+<j>)    (3-4) 


where  <j>  -  Cos"   <g 

The  oscillatory  nature  of  the  C(t)  curve  in  Fig.  10 
is  of  this  form.  When  the  control  valve  metering  area  is 
small  (say  A^  =  .0194  in2  as  in  Fig.  2)  P_(t)  deviates 
considerably  from  being  a  step  input.   The  result  is  that  the 
C(t)  curve  shows  no  overshoot  at  all. 


E.   Rise  time  Consideration  for  Solenoid 
Control  Valve 


One  more  thing  should  be  mentioned  here.   The  3-way 
solenoid  on-off  type  control  valve  discussed  in  this  report 
was  assumed  to  open  instantaneously  at  time  t  =  0+.  In  the 
actual  case,  a  finite  rise  time  is  required  for  the  solenoid 
valve  to  fully  open. 

area=A0  -^  area=A0 


S 


0 


Assumed 


rise  time  ■  0 


~7 


o 


"H  rise  time 


Skinner  solenoid  control 
valve  rise  time  =  0.0045 
sec. 


Fig.  11.   Rise  Time  of  Solenoid  Control  Valve. 
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This  rise  time  will  vary  considerably  according  to  the  fluid 

media,  and  the  applied  voltage.  As  a  typical  example,  a 

Skinner  Electric  Co.  V-5  series  three-way  solenoid  valve 

requires  .0045  sec  (AC  operation)  to  fully  open  if  the  fluid 

is  air.   The  response  times  calculated  in  this  report  would 

probably  be  .003  to  .004  sec.  greater  in  the  case  of  an 

actual  solenoid  valve  and  actuator.   Even  so,  they  would  still 

be  0.003  sec.  less  than  the  0.033  sec.  response  time  desired 

for  2  cycle  breaker  interruption. 

• 
F.   Selection  of  At  for  Computing  Response 
Characteristics 

In  numerical  approach  to  the  solution  of  the  response 
characteristics  -of  a  nonlinear  system,  the  proper  choice  of 
the  time  increment  At  is  very  important.   The  size  of  At 
for  computation  is  determined  by  the  system  parameters  and 
the  degree  of  accuracy  demanded.   Choice  of  a  large  (coarse) 
value  for  At  will  cause  large  error  in  the  computated  results, 
while  choice  of  a  small  At  will  require  more  computing 
iterations.   Initially,   t  =  0.0001  sec.  can  be  selected  for 
a  first  trial.   If  the  response  curve  for  Pc(t)  comes  out 
with  peculiar  oscillations  of  relatively  high  frequency,  this 
indicates  that  choice  of  a  smaller  value  for  At  is  required. 
Usually,  this  kind  of  oscillation  only  occurs  over  certain 
time  intervals.   Therefore,  the  computer  program  was  written 
such  that  the  smaller  value  of  At  was  only  used  in  the 
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effected  portion  so  as  to  save  unnecessary  computing  itera- 
tions.  For  instance,  in  computing  the  results  of  Fig.  10, 
a  small  value  for  At  (0.00001  sec.)  had  to  be  used  in  the 
time  interval  0  to  0.032  sec. 
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CHAPTER  IV 
CONCLUSIONS 

In  this  report,  the  response  characteristics  of  a 
nonlinear  compressible  fluid  actuating  system  were  investi- 
gated.  Fluid  actuating  system  of  this  type  are  used  for 
controlling  the  operation  of  air  blast  circuit  breakers  used 
to  protect  EHV  power  transmission  systems.  Determining  the 
response  characteristics  of  such  a  system  is  quite  involved 
because  they  are  nonlinear.   Use  of  approximate  techniques 
to  investigate  the  response  resulted  in  the  following 
conclusions: 

1.  Numerical  techniques  can  be  employed  to  obtain 
the  approximate  response  characteristics  of  such 
systems.   The  accuracy  of  this  approach,  approximately 
2%,    is  acceptable  for  engineering  design  purposes. 

2.  The  analysis  made  in  this  report  is  based  on  certain 
assumptions,  therefore,  some  deviation  in  the  per- 
formance will  be  expected  in  actual  case.'  Specifi- 
cally, the  neglecting  of  the  rise  time  of  the  solenoid 
valve  may  cause  some  errors. 

3.  The  response  time  of  the  actuator  is  strongly  in- 
fluenced by  the  size  of  metering  area  of  the  control 
valve. 

4.  The  supply  tank  pressure  drop  occuring  during  a 
complete  opening  operation  did  not  appreciably 


33 


effect  the  response  time  for  up  to  3 .in.  of 
actuator  travel,  for  supply  pressure  drops  as  large 
as  36  psi. 

5.  A  large  control  valve  metering  area  caused  Pc(t) 
to  become  approximately  a  step  input.   This  re- 
sulted in  overshoot  and  oscillation  in  the  C(t) 
response  curves. 

6.  -A  limit  for  increasing  A0  in  an  attempt  to  make  tr 

small  exists.   In  the  system  discussed  in  this 
report  when  A0  was  increased  from  0.25  to  1.0  the 
decrease  in  tr  was  only  0.0037  sec. 

7.  Proper  choice  of  the  time  increment  &t   used  for 
computing  the  results  is  important.   If  a  high 
frequency  oscillation  appears  in  the  computed  re- 
sults for  Pc'(t),  a  smaller  At  should  be  tried. 

8.  The  system  (with  Vt  =  1000  in3,  AQ  =  0.0194  in2) 
presented  in  this  report  has  a  response  time  of 
0.0225  sec.  for  3  in.  of  piston  travel.   If  0.004 
sec.  is  allowed  for  the  rise  time  of  an  actual 
solenoid  valve,  the  response  time  0.0265  sec.  would 
still  satisfy  the  interrupting  time  of  2  cycles. 
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APPENDIX  A 
FLOW  CHART  OF  COMPUTER  PROGRAM. 


36 


Eh 
+ 

E-i 


N=  2 


N=N+1 


K=l 


COMPUTE 
CD(lO 


YES 


O 

o 
o 

o 

+ 

Eh 


i 

r 

COMPUTE  PTD  PCD  CDD 

1 

K=K+1 

NO 


37 


APPENDIX  B 
FORTHRAN  PROGRAM 


NONLINEAR  ANALYSIS   OF  THREE-WAY  VALVE  AND   CYLINDER  OPERATION       ■  ' 
S?Ss^N  PT(lOOO),PC(lOOO),CD(lOOO)>C(lOOO),RT(lOOO)>RC(100O) 

DIMENSION  W(IOOO)        ■ 

DIMENSION  PTD ( 1000 ), PCD ( 1000 ),CDD( 1000) 

200  ™T(1^20^^ 

READ  (1,200)  PCD(l5,CDD(l),CD(l),C(l),RT(l),RC(l),C0 

95   FoSl/5X^ 
WRITE   (3,95) 
T=0. 00001 
DO  10  N=2,1000 

Wtt-1 

PTD(M)=PTD(M) 
PCD(N)=PCD(M) 
CDD(N)-CDD(M) 
CD(N)=CD(M) 

K»l 

PT(N)=PT(M)  +  (PTD(N)+PTD(M)*DT/2. 

PC(N)=PC(M)+(PCD(N)+PCD(M)*DT/2.     .  ■    • 

CD(N)=CD(M)+(CDD(N)+CDD(M)*DT/2. 

IF(N-5)25, 25, 88 

IF(CD(N))131,131,25       ■      ■ 
C(N)=C(M)+(CD(N)+CD(M)*DT/2. 

RT(N)=RT(M)*(PT(N)/PT(M),,**0.714 

RC  ( N  )=RC  ( M )  *  ( PC  ( N ) /PC  (M ) **0 .  714 

CR=PC(N)/PT(N) 

IF(CR-0. 523)30,30, 32 

W ( N ) =DK*A  0*SQRT ( 180 . 5 *RT ( N ) *PT ( N ) ) 

GO  TO  60 

S=1.-(PC(N)/PT(N))**0.286 

IF(S)40,33,33 
W(N)=SQRT((S)*((PC(N)/PT(N))**1.428)) 

W( N )=51 . 9*DK*A0*SQRT (RT (N ) *PT ( N ) ) *W( N) 

GO  TO  60 

S=l . - ( PT (N ) /PC ( N ) ) **0 .286 

W(N)--SQRT((S)*((PT(N)/PC(N))**1.42d)) 

W(N)=51.9*DK*AO(SQRT(RC(N)*PC(N))*W(N) 

PTD(N)=-1.4*PT(N)*W(N)/(TV*RT(N))  ,    %11 

PCD(N)=1.4*(V/(N)-RC(N)*A*CD(N))*PC(N)/(A*RC.(N)*(C0+C(N))) 

CDD(N)=386.0*((PC(N)-PC(l))*A-BP>^D(N)-SK*C(N))/(PM+CM) 

K=K+1 

lF(K-3)    20,20,70 
70  WRITE    (3-80)   T,PT(N),PC(N),CD(N),C(N) 

IF(N-650)10,10,12 
12  DT=.0001  ' 
10  T=T+DT 
131  STOP 
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This  report  presents  a  nonlinear  analysis  of  a  three- 
way  control  valve  and  cylinder  operating  with  a  compressible 
fluid.  The  system  consists  of  a  supply  tank,  a  piston  type 
actuator,  and  an  on-off  type  solenoid  control  valve.   Such 
actuating  systems  are  used  to  operate  air  blast  type  circuit 
breakers  which  are  used  to  protect  power  transmission  net- 
works. 

The  characteristics  of  compressible  fluid  flow 
through  valves,  plus  displacement,  velocity,  and  accelera- 
tion equations  are  derived..   Then  equations,  all  nonlinear 
except  one,  are  solved  by  numerical  method. 

Response  curves  for  actuator  displacement,  actuator 
velocity,  and  supply  tank  pressure  at  functions  of  time 
were  calculated  and  are  presented. 

The  results  indicates  that  the  response  time  of  the 
actuator  is  most  strongly  influenced  by  the  size  of  the 
metering  area  of  the  control  valve.   It  was  found  that 
supply  tank  pressure  drops  up  to  36  psi  did  not  appreciably 
effect  the  response  time  for  actuator  travels  up  to  3  inches, 

The  results  also  indicate  that  a  large  control  valve 
metering  area  will  result  in  overshoot  in  the  response 
curves  if  the  damping  ratio  for  the  actuator  is  less  than 
1.0. 


